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Abstract 

Using AdS/CFT, we calculate the energy density of a quark moving through a thermal 
state of A/" = 4 super- Yang-Mills theory. Relying on previous work for momentum-space 
representations as well as asymptotic behaviors, we Fourier transform to position space and 
exhibit a sonic boom at a speed larger than the speed of sound. Nontrivial structure is found 
at small length scales, confirming earlier analytical work by the authors. 
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1 Introduction 

Heavy ion collision experiments at the Relativistic Heavy Ion Collider (RHIC) are believed 
to probe properties of a deconfined state of hadronic matter, the Quark Gluon Plasma 
(QGP) [H El El H] . The phenomena of jet-quenching and jet-splitting observed at RHIC are 
challenging to understand theoretically. A review of theoretical treatments of jet-quenching 
can be found, for example, in [5]. Theoretical works on jet-splitting include [HI [71 [H El |TOl 
[TT] . Clearly it is desirable to study similar phenomena in strongly coupled A/" = 4 super- 
Yang- Mills theory (SYM), where substantial progress can be made via AdS/CFT starting 
from calculations in supergravity. While direct comparisons between SYM and QCD are 
fraught with uncertainties, such calculations at least provide an independent perspective 
on dissipative effects arising from a quark propagating through a strongly coupled thermal 
medium. In this paper the moving quark is infinitely massive, pointlike, and fundamentally 
charged. 

A heavy quark moving at constant velocity in an infinite, static, thermal bath of A/" = 4 
super- Yang-Mills theory can be described in the dual gravity picture by a string whose 
endpoint is on the boundary of the AdSs-Schwarzschild spacetime. This string is responsible 
for the drag force on the quark [121 [13] (see also the closely related work [T3|, as well as other 
work [15] on jet-quenching in AdS/CFT). The configuration mentioned has a steady-state 
approximation built in: the trailing string is moving at the same velocity v as the quark. 
Our convention is to take this velocity to be in the Xi direction: v = (f, 0,0). We employ 
mostly plus metric, so that = Xi but x^ = —Xq = t. 

The string perturbs the geometry of AdSs-Schwarzschild in a manner explained in |16j . 
From the asymptotics of those metric perturbations one may extract the expectation value 
of the gauge theory stress-energy tensor. In [TH] the Fourier coefficients of this stress-energy 
tensor have been evaluated. It is clearly desirable to pass to a position space description: this 
is a necessary step before hadronization can be considered, and it also helps one's intuitions 
about the dissipative mechanisms involved. The main result of this paper is the computation 
of the energy density at a fixed time as a function of position x. More precisely, we compute 
a dimensionless quantity 



Here (Too(0, x)) is the thermal expectation value of the total energy density in the gauge 
theory at time t = 0. For convenience, an overall factor of yl — v^ has been brought out 




(1) 
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explicitly in ([T]). We always use spacetime coordinates such that the plasma is at rest. The 
second term inside the parentheses in ([T]) is the contribution to (Tqq) from the thermal bath. 
We will generally work with dimensionless position space and Fourier space variables: 



X = ttTx, K = k/nT . (2) 

With the Fourier coefficients £{K) in hand, the problem we face is simply to compute a 
three-dimensional Fourier transform: 

A difficulty in carrying out this Fourier transform is that S{K) grows linearly with K and 
is also singular at small K. Our strategy will be to find asymptotic expressions with known 
analytic Fourier transforms at large and small wave-numbers, subtract them away, and pass 
the remainder through a fast Fourier transform (FFT). The simplest of these expressions is 
the Coulombic near-field of the quark, which is proportional to in the rest frame of the 
quark and takes the form 



(1 - v^f/^ [Xf + (1 + v^)Xj 
127r2 [Xf + (1 - v^)Xl 



C-CoulomblAj - _ ^ , ON ..2l3 



in the rest frame of the plasma. Here X\_ = a/X| + X|. Finite-temperature corrections to 
(jlj) which are still singular at X = were found in [T71 [T8l [19] and are developed further 
in this work. These corrections may be interesting in their own right, because they provide 
some gauge-invariant information about the energy loss mechanisms close to the quark. They 
might be used to seed hydrodynamical simulations. 

In presenting our results we will usually consider the quantity 

E{X)=S{X)-ScouioMX) (5) 

/ total energy density\ energy density 

= I energy density - of plasma ] - of moving quark . 
V of system in equilibrium/ in vacuum 

Up to an overall prefactor (see ([1])), E{X) is the total energy density of the moving quark 
and the thermal bath, minus the energy density of the thermal bath in the absence of the 
moving quark, minus the energy density of the moving quark in the absence of the thermal 
bath, all at a fixed time. We will refer to E{X) somewhat loosely as the net energy density. 
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It may be positive or negative. 

The reader wishing to skip over technical details can find a summary of our results in 
section 12.31 Briefiy, we find good agreement both with the analytical estimates [TTl [181 EH] 
at small length scales, and with linearized hydrodynamics at large length scales. 

During the late stages of this project, we learned of a similar study of the energy density 
which has some overlap with this work |2Uj . 



2 Real space calculation of the energy density 



In earlier work [161 lIH]! the relation between the thermal expectation value of the stress- 
energy tensor T^^^, and the metric fiuctuations due to the trailing string has been worked 
out in detail. Here we present a brief summary. Let us write the spacetime metric as 
ds'^ = G ^ydx^dx" , with G^^y = Gf} + h^^y, where G^S is the AdSs-Schwarzschild metric 



^4) = G^S^d^'d^"" = ^ 



-g{z)dt^ + dxl + g{z)-^dz^ 



9{z) = 1 - 



4 ' 




and h^v is the first order correction to it in response to the trailing string. The temperature 
T of the field theory is related to the location zq of the black hole horizon through zq = 1/ tcT. 
The stress-energy tensor of the boundary theory is proportional to the fourth order coefficient 
in the expansion of the metric fiuctuations h^y{z) at small z. The stress-energy tensor is 
traceless, as expected for a conformal theory, but it is not conserved. The non-conservation 
simply means that energy, as well as momentum in the direction of the quark's motion, are 
deposited into the thermal plasma at a constant rate [16]. This rate reproduces the drag 
force as computed in [121 [13] • 

In order to obtain the stress-energy tensor explicitly, one has to solve the linearized 
Einstein equations for the metric fiuctuations h^p sourced by the trailing string. These 
equations take the form 

V'^P'^hp, = J""" , (6) 

where V^'^^^ is a second order differential operator (a Lichnerowicz operator) and J'^'^ is the 
source term generated by the string. These equations can be reduced to five decoupled second 
order differential equations written in terms of gauge-invariant quantities called "master 
fields" [IH]. Because of parity considerations, two of these equations have no source terms. 
The two corresponding master fields can be set to zero, and so they don't contribute to hfj,iy. 
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The remaining three master equations can be written in the form 



^d,z-^gd-, - 1 + v^kl {g-' - 1) + ~eVx{z~k)\ ~k-'^x{z~k) = Jx{z~k) , (7) 

with X = S, V, and T corresponding to "scalar," "vector," and "tensor" perturbations, 
respectively (all with even parity). Also, (fci, /c2, k^j are the variables conjugate to {xi, X2, x^) 
(recall that our metric signature is mostly plus), and k"^ = k"^ — v'^kf. In what follows we 
shall find it convenient to switch to the dimensionless variables K = z^k and X = x/zq, as 
defined in ([2]). We will also use a convenient rescaling of the "depth" coordinate in AdSs, 
namely Z = zk. 

The quantity S{K) appearing in (j3]) can be obtained from a small Z expansion of the 
scalar master field $5: 



$s = . . . - R-"^ (^-12nS{K)vKi + z6i;^) Z^ + 



where K = y/Kf{l - v^) + Kl and = ^Kl + Kl (see [IE] for details.) 

The small K asymptotics of E,{K^ may be obtained by expanding $5 in power series in 
J^", $5 = ^ K^ipni and solving the corresponding equation in ([7]) perturbatively in small K 
[TB| [TB] (see equation ([TB|) below.) The large K asymptotics can be obtained in a similar 
manner as a power series in 1/ K, which may also be regarded as a series in positive powers 
of T. This series was calculated up to order 0{K~^) in [191 118]. In appendix \K\ we rederive 
these results and extend them to order 0{K~^). The values of S{K) for intermediate values 
of K have to be evaluated numerically. We do this using the method developed in [TB] . 

With small K and large K asymptotics in hand, we consider a decomposition 

£{K) = £ijy{K, /iuv) + £ir{K, /iir) + £res{K, /iuv, /im) • (9) 

Here 8\jy{K, fiuy) agrees with the large K asymptotics up to and including terms of order 
K~^; SiniKjfJ'm) almost agrees with the small K asymptotics up to and including terms of 
order (see equation (12T!) and the discussion following it for more precise details); and 
SresiK , fixjy , fii^) is uniformly bounded and integrable. The parameters /iuv and /iir are 
pure numbers which can be adjusted to make the residual part Sj-csiK) as small as possible. 
The precise analytic forms of S\jy{K , fiuy) and SiR{K,fLm) will be explained in sections ETT] 
and 12.21 respectively. 

Ideally, £\jy{K, fiuy) and £^iR(i?, /^ir) should admit analytic Fourier transforms. We 
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didn't quite realize this goal: the Fourier transform of £i^{K, /iir) is left in the form of a one- 
dimensional Fourier integral which must be performed numerically. The Fourier transform 
of S^cs{K, fi\jy, /iir) must also be performed numerically, via a three-dimensional fast Fourier 
transform. In section 12.31 we report the results of numerics for three values of velocity: 
V = 0.25, V = 0.58, and v = 0.75. The second of these is only slightly larger than the speed 
of sound in the thermal plasma, Cg = 1/ Vs ~ 0.577. 

2.1 Near-field asymptotics 

Using methods explained in section [2] and appendix El one finds 

S{K) = sSiiK) + 0{K-') , (10) 

where 



'in 



24:K'' 



where, as before, k = ^Kf{l -v^) + Kl and K± = ^K^ + Kl The difference S{k) - 
Sljy{k) is small at large K but large at small K because of the inverse powers of k that 
appear in (fTTl) . This is bad because our eventual aim is to find an analytic approximation 
to S{K) that is good both in the UV and the IR. The bad IR behavior of ( fTTl) can be cured 
by using the identity 



-n/2 

(12) 



and expanding the quantity in square brackets to just enough terms to keep the 0{K ^) 
accuracy that 4v(^) 

possesses in the first place. Through this procedure one obtains 



tuY[K,^UY) = — \ , ^ \- ■■■ , [IS) 

^4 48 + 

where the omitted terms are polynomials in Ki, /iuv, and v times negative powers of 



K"^ + /iuv- These terms are straightforward to work out, but their precise form is long 
and not very enlightening. 
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The Fourier transform of Suy{K, fiuy), as well as several other Fourier transforms re- 
quired in later sections, can be worked out starting from 



and taking appropriate derivatives of it. Here Kj,(z) is a modified Bessel function of the 
second kind. For the case at hand we use f|T4l) with c? = 3 and K ^ K. 
The Fourier transform of £-[j\{K, /luv) is 

2 

£uv(^, /iuv) = ^""^ 3 - f/"uvX (2v^Xl + (-2 +v^ + v^) X^) Ki(/iuvX) 

967r2(l -^2)2X4 V V ' ' ^ 

+ 4 (2v^Xl + (-1 + i;^) X^) K2(^uvX)) + . . . , (15) 

where we have defined X = ^JXl/{l — u^) + X\. The terms written explicitly in (1151) 
correspond to those appearing explicitly in f|T3|) . Additional terms have similar forms, and 
it would not be very illuminating to write them out explicitly here. 



2.2 Far-field asymptotics 

As described at the beginning of this section, the IR asymptotics of the energy density may 
be obtained by solving equation ([7]) perturbatively in small K. One obtains 

- ?AK,v (v^ + 1) 3A-;«^ (^1(2 + v') + 2Kf(l + v')) 
^<''' = -2.(A-(l-3.^) + Al)+ 2Mm-3v') + Kl)' + "'''^ ' ''"^ 

A problem with f|T6l) is that the Fourier transform of negative powers of Kf{l — Sv"^) + Kj_ is 
not well defined when v"^ > 1/3 due to the poles which appear on the real axis. As suggested 
in [16], the quantity 

^(resummcd) / tTJ-n _ 

^IR — , rr^ n2 U'J 

Kj^ + m(Ai) 

with 



271 {1 + v"^ - tKiv{2 + v^)) 
{3v^ - 1)(1 + v^)Kf + 2iKfv{l - iP) 



1 + - iKiv{2 + v^) 

is a more uniform approximation to S{K) for small -fC than the terms shown explicitly in ( fT6l) . 
The expressions (fT^ and (IT7I) agree up to C(-ft') corrections, but the denominator of (fTTI) 
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shifts the poles away from the real axis. We describe this modification of the denominator as 
a "resummation" of the series (1161) because a single rational expression, fll7p . includes both 
terms in f|T6|) FI 

The problem now is that an analytic Fourier transform of (^^{^ Yeal space is 

unavailable (as far as we know) because of the cubic terms in the denominator. However, 
the integration over K± = {K2, K3) may be done analytically using (ITll) with d = 2 and 
n = 1. The result is 

sSr-'^'m = / ^^e'" ^<^';' = r !^^e"<""^K,(m(K,)X,) (19) 



where m{Ki) = \Jrn{Kif is chosen so that Re{m(A'i)} > 0. The last integral in (1191) 
does not seem to be expressible in closed form. However, for X±_ > 0, the integrand has 
exponential decay for large \Ki\, and it is continuous. So it can be handled reliably with 
numerical methods, i.e. a one-dimensional FFT. 

The large K behavior of (^^{^ jg (9(1//^^), so it brings back problems with 

the UV region that £^uv(-^; A^uv) was supposed to have fixed. Let us therefore consider the 
following corrected form: 

The Fourier transform of the second term in (l20i) is similar to ( IT9i) . At small K one finds 

S{K) - S,n{K, /xir) = + 0{K) . (21) 

Better would have been a form for SiYi{K, fim) which matched S{K) up to 0{K) corrections, 
with no (9(1'^°) term as in (!2T!) : for instance, some exponential suppression of this term could 
have been arranged. But, because the deviation of the subtracted quantity Si^{K, /xir) from 
S{K) is finite at small K, it does not introduce divergences at small values of the momenta 
and its Fourier transform can be easily dealt with numerically. 



^As will be discussed in section [31 this resummation can be motivated physically in terms of passing from 
inviscid to viscous hydrodynamics. 
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2.3 Numerical results 



We obtained S{K) numerically using minor improvements over the code used for [TB]. We 



chose yUuv = fJ'iR = 1- A three-dimensional FFT of £^rcs = ^ — ^vy — "was performed 
on a grid with 128 points on each side and with each component Ki running approximately 
from —10 to 10. One-dimensional FFT's of Sj^i were carried out on grids with approximately 
3000 points, with Ki again running approximately from —10 to 10. Edge effects can be 
expected in position space near the boundary of an FFT. To smooth them out we multiplied 
the position space results of FFT's by a approximation of the top-hat function, i.e. a 
function which is outside a large cube, 1 inside a smaller cube, and has one continuous 
derivative everywhere. 

In figure [1] we plot the net energy density E = £ — f^couiomb as a function of Xi at 
various velocities and various values of X±. To get a feel for the scales involved, note that 
if T = ^ GeV ~ 318 MeV, then X = 1 corresponds to a distance 1 GeV~^ ~ 0.2 fm from the 
quark. 

The lower left plot in figure [1] shows that our infrared subtractions are not as accurate 
multiplicatively for f < 1/ as they are for > 1/ ^/S. Probably a better subtraction 
scheme could be contrived for v < starting from the quantity £^(^^^*^"<^) discussed in 

section [3l Note however that the magnitude of E is very small in the region in question. The 
middle left plot in figure [1] shows that the asymptotic form Suy{X , fiuy) deviates strongly 
from the full expression for E on intermediate length scales when v = 0.25. This is not 
particularly alarming: the scales in question are far from the ultraviolet regime. However, 
one might hope by cleverer choice of asymptotic forms (perhaps a more felicitous choice of 
yU-uv and /iipt) to obtain a more uniformly accurate approximation to E{X). 

For supersonic velocities, one can compare the structure associated with a Mach cone 
to the scale of hydrodynamical broadening. Ax = y^Tgt. Here Tg is the diffusion length for 
sound waves. We plug in t = ^csc6'jv/, where 9m = arccos ^ is the Mach angle, because 
this is the time it takes for sound to travel from the quark to a point on the Mach cone a 
transverse distance x± from the axis of the quark's motion: see figure [21 The sound waves 
propagate at the Mach angle, whereas our plots in figured] are along the Xi direction. Thus 
the expected broadening in xi is 




(22) 
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Figure 1: (Color online.) E = S — i^couiomb as a function of dimensionless position Xi for 
various values of offset Xp = X± from the quark and velocity v of the quark. In each plot, 
the black curve is E. The red curve is The green curve is S^es with smoothing (see 
the main text). The blue curve is £\jy — £^couiomb- The dotted purple curve is from inviscid 
linearized hydrodynamics (see the discussion in section [3]). The green dot corresponds to 
the Mach angle. The blue dots show a width typical of hydrodynamical broadening (see the 
main text). The red dots show the points where the smoothing of S^-cs starts. The quark is 
at Xi = X^ = 0. 
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Figure 2: (Color online.) An estimate of hydrodynamical broadening of the Mach cone. The 
gray line, offset by a distance from the quark's trajectory, is the axis along which 

we plot E = £ — £'couiomb in figure [1] The green dot in that figure is at the intersection of 
the gray and green lines in this one; likewise the blue dots are at intersections of the gray 
line with the blue lines. 



Using Cs = l/"\/3, = l/SvrT, and Xi = vrTxi, one finds 



AXi = ttTAxi 



(23) 



Figures [31 HI and \5\ show contour plots of the net energy density E = S — £^couiomb 
for various values of the quark velocity. The structures exhibited in E = S — i^couiomb are 
notably scale-dependent. Observe for example from figure[3]that at f = 0.25, E > for small 
negative and large positive Xi and E < for small positive and large negative Xi. Recall 
that if T = ^ GeV ~ 318 MeV, then X = 1 corresponds to a distance 1 GeV~^ ~ 0.2 fm from 
the quark. 

Note that the double-peaked feature in the upper right plot in figure [1] corresponds to 
the structure near the origin of the top plot in figure [5l As observed in [T8l [T9] , in three 
dimensions this structure is a forward lobe and a backward-leaning cone of energy over- 
density (as usual, relative to the sum of the energy density of the moving quark and the 
thermal bath in the absence of interaction between the two), with regions of under-density 
in complementary regions. 

In the large scale plots one sees that the net energy density falls quickly to zero for 
positive Xi when v > 1/ We suspect that £ (not E) falls exponentially in this direction. 
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E for v=0.25 




-0.05 0.05 



E for v=0.25 




Xi 

E for v=0 .25 




Figure 3: (Color online.) Contour plots oi E = £ — ^^couiomb as a function of dimensionless 
position coordinates Xi and Xp = X±, for v = 0.25. Note that the top plot shows the smallest 
structures while the bottom plot shows the largest. Orange and red regions correspond to 
E > 0; white regions correspond to E' ~ 0; and blue regions correspond to < 0. The 
energy density of the thermal bath is not included in E. The three-dimensional energy 
density profile is axially symmetric around the Xi axis. The black dot is the position of the 
quark: Xi = X± = 0. 
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E for v=0.58 




E for v=0.58 



E for v=0 . 5f 




Figure 4: (Color online.) Contour plots oi E = S — £^couiomb as a function of dimensionless 
position coordinates Xi and Xp = X±, for v = 0.58. Note that the top plot shows the smallest 
structures while the bottom plot shows the largest. Orange and red regions correspond to 
E > 0; white regions correspond to E' ~ 0; and blue regions correspond to < 0. The 
energy density of the thermal bath is not included in E. The three-dimensional energy 
density profile is axially symmetric around the Xi axis. The black dot is the position of the 
quark: Xi = X± = 0. The dashed green line shows the Mach cone. 
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E for v=0.75 
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Figure 5: (Color online.) Contour plots oi E = £ — £^couiomb as a function of dimensionless 
position coordinates Xi and Xp = X±, for v = 0.75. Note that the top plot shows the smallest 
structures while the bottom plot shows the largest. Orange and red regions correspond to 
E > 0; white regions correspond to E' ~ 0; and blue regions correspond to < 0. The 
energy density of the thermal bath is not included in E. The three-dimensional energy 
density profile is axially symmetric around the Xi axis. The black dot is the position of the 
quark: Xi = X± = 0. The dashed green line shows the Mach cone. 
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For negative Xi and fixed X±, E decays as 1/Xl (see equation (128|) ). 



3 Conclusions 

Through a combination of analytic and numerical methods, we have calculated the position 
space energy density of a quark moving through a thermal state of A/" = 4 plasma. More 
precisely, we considered the net energy density, obtained by subtracting away the constant 
contribution from the thermal bath as well as the field the quark would have generated in 
the absence of the bath. This net energy density provides some gauge-invariant information 
about energy loss. It has interesting structure on multiple scales, from sizes ~ 0.01 fm to 
~ 2 fm if we take T = 318 MeV. By comparing asymptotic forms (analytic or semi-analytic) 
to numerical results we can be reasonably confident that we have accurately characterized 
these structures. The power of AdS / CFT is that all length scales can be treated at once. 

The main features of our results can be understood in terms of analytic approximations. 
At large distances and supersonic velocities, the Mach cone can be described by the small 
momentum asymptotics of the energy density, 

as found in [16]. (We have excluded from the denominator in flM|) some cubic terms that 
vanish on the Mach cone.) If the cubic term in ^IM were dropped, then for v > l/VS 
there would be poles on the Ki axis corresponding to zeroes of — Sv^Kf. An integration 
prescription is required to pass these poles, and the physical one is given by the —ivK'^Ki 
term. Carrying out this contour integral one finds a Mach cone structure in real space, as 
could have been expected from linearized hydrodynamics. 

The relationship of fl2^ to linearized hydrodynamics is worth understanding a bit better. 
According to linearized hydrodynamics, the energy density obeys an equation of the form 



^2 Q2 

^-^(c' + r.^t) e = sources. (25) 



(See for example [H].) Comparing the Fourier transform of (l25l) to the denominator of 
( |2^ . and recalling that all dependence of e on time and xi is in terms of xi — vt, one 
finds perfect agreement if = 1/3 and F^ = l/37rT. In general, F^ = Arj/SsT, so this 
provides an independent check of the result rj/s = 1/4tt [21]. Evidently, the cubic term in 
the denominator of fl2^ arises from viscosity. If it is removed, corresponding to sending 
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?7 ^ 0, then the energy density becomes singular on the Mach cone. More exphcitly: the 
inviscid hnearized hydrodynamics approximation to fll7p is 

C (inviscid).. 7^ ^ivjl + V^) Ki 

where the infinitesimal positive quantity e provides the correct pole passing prescription. 
One easily finds 

c, (inviscid), ^^ _ 3f (1 + V^) Xi 

^ '~ 87r2 (X2 + (1 - 3t;2)Xi)3/2 

ioT V < 1 / y/S, and 

,,™..,,^.) ^ ^"'^7^',^..^^^:\,,^^.^3/. ^(-A. - X,VW-1) (28) 

for V > 1/ y/S. The functions fl27j) and fl28|) are shown in the large scale plots in figure [TJ 
The result ( !28|) can be supplemented by a singular distribution of positive energy right on 
the Mach cone, but it is best regarded as an asymptotic form for large negative Xi. 

The charm of the AdS / CFT calculation is that it interpolates smoothly all the way from 
down to the Coulombic near-field of the quark, for which e ~ in the quark's rest 
frame. A good understanding of the structure observed at scales x ^ 1/T can already 
be obtained from considering just the analytic (!?(T^) correction to the Coulombic field 
p!8| [T9] . In particular, the over-density of energy ahead of the quark for v = 0.25 and 
V = 0.58 and the more complicated structure near the quark for v = 0.75 were predicted 
by that work. These non-trivial structures illustrate the need (already well recognized) to 
supplement hydrodynamics with a physical prescription of how energy is deposited near the 
quark. The next order O(T^) corrections have been calculated in appendix Rl 

Extensions of the analytical computations in appendix |Al as well as numerical treatments 
of other components of the stress tensor, are underway. When complete, perhaps they will 
provide an interesting alternative perspective to perturbative intuitions about energy loss 
mechanisms in QCD. Caution is appropriate when comparing AdS/CFT calculations to 
QCD in this context because A/" = 4 SYM is conformal, and therefore equally strongly 
coupled at the smallest and largest of scales. 
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A Large momentum asymptotics 



The starting point for computing the large momentum asymptotics of the energy density 
£{K) is the scalar master equation in ([7]), with the potential 



12 {2Z^ (l + a2)+2ir4 {l + a^) + Z^ {-12 + K\l 
(6Z2 + A'4 (l + a2))2 



a 



and source term 



(1 + a^){- 54:Z^a + 3K^a (l + a^) 



a 



- 6iZ^K^ (2 + lla^) + 3Z^K^a (5 + lla^) + SiZ^K^ (2 + Sa^ + 3a^) 
+ Z^K^a (90 + /s:^ (-2 + a^) (l + a^)) 



{ - 162tZ^K^a^ + {9t) Z^K^a^ (l + a^) + 3A^^ (-2 + a^) (l + a^) 
- 18Z^ (-2 + 7a^) + 3Z^A'^ (2 + Ta^ + 23a^) 

+ Z'K^{-2 + a') (90 + A'^(-2 + a2) (1 + ^2)) 



(29) 



(30) 



as constructed in [181. Here 



Ai 



, AT-iZ , K + Z 
log ~, + « log 



K + iZ 



K -Z 



(31) 



and a = vKi/K. Equation (I7j) can be solved perturbatively at large K using methods 
developed in [I71[l8]@ 

Writing $5 = Y2K~'^ipn and solving the equations of motion perturbatively for each tpn, 
one finds 

^0 = l^rZ^C, + C2)(Lo(Z) - Io(Z)) + C2Z , (32) 
where Iq is a modified Bessel function of the first kind, Lq is a modified Struve function, and 



^In [T7] the large momentum asymptotics were computed to order 0{T'^) using the WKB method which 
required a piecewise approximation to the Schrodinger potentiaL This approximation breaks down at order 
©(T^) so instead we use an iterative scheme developed in [l^ for extracting the corresponding correction. 
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Ci_2 are given by 

a (-2 + a2) (2-a2 + t;2 (1 + a^)) _ 3a (2-a2 + y2 (1 + ^2)) 

1 + 1 + 

At the next order one has 

^2 = C^Z^ + , (34) 

with 

C4 = -izCitt C3 = (-5 - 2 + (1 + (9 + 2 a^)) . (35) 
'?/'4 can be found using the Green's functions method. Consider 

- = Z^Ko{Z) / dxx"%{x)J{x) + Z%{Z) / dxx-^Koix)J{x) , (36) 
Jo 



where now 



and 



J{x) = -27r(Ci + C2)Z'A,{Z) + Ao{Z) ^ 5,Z' + e.^' (37) 



12 (Ci + C2) 71 



2 (1 + ^2) 

6 (-4C2 (l + a2) + 9a (2 - ^2 + ^,2 (1 + ^2))) 



(38) 



1 + a2 

64 = A (Ci + C2) Tx (39) 
, (C1 + C2) vr (-l+a^) 



(40) 



= ^ '-^ ' ,^ ^ \.2 ^ (41) 

(1 + a^) 

-6 ('2C2(l + a2)2_ct (8- lla2- 10a4 + t;2 (7 + 11 ^2 + 4 ^4))) 

es = ^ 2 ^ (42) 

(l + a2)2 V ; 

_ a (60 -8^2 ^7^4^ ^2 (^g^^^ ^2 _ 7^4)) 

6(l + «2) ^^2) 

^ -{o^' (-2 + c.^) (2-c.2 + ^;2 (1 + ^2))) 
' 18(l + a2) • ^ ^ 

Since for small values of their arguments J ~ x^, Iq = 1 + O (^^) and Kq ~ 1 + InZ + 
C(Z2), one has 

_ (-V2^ ^ ^1^^ ^ (-^^(C^i + C^2) + C3 - i rfxx-^Ko(x) J(x)) + 0{k-') 

(45) 
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where the upper hmit of the second integral in (l36l) has been taken to infinity. 
Carrying out the integral, one finds 



^ ttK (2 - + (1 + a^)) a (-5 - 2 + (11 + 2a^)) 

-nt = — h I = 

24 ISK 

vr(7 + 3a^) {2 - + {1 + a^)) ^ ~ 
24^3 

B Numerical Fourier transforms 

In this appendix we describe our conventions for numerical Fourier transforms and how they 
relate to the underlying code, namely Mathematica's implementation of the FFT. Although 
the contents of this section are entirely elementary, we found it a useful reference. Using 
the default settings, Mathematica's FFT code transforms a list u^. of complex numbers with 
1 < r < into another list according to 



^^M,e2-(^^-^)(^-i)/^. (47) 



r=l 



(We discuss here only the case of a one- dimensional FFT; higher-dimensional cases follow in 
the expected fashion.) Consider now a real function 



°° '^^ ikx . 



m = / ^e-V(A^) (48) 

which is peaked at x = xq and decreases to on a length scale Ax. The Fourier transform 
f{k) satisfies f{—k) = f{k)*, and for smooth /(x), it decreases to roughly on a scale 
Ak = 1/Ax. Taking this as a defining relation for Ak, let us further define 

Sk=J^Ak 6x=J^Ax. (49) 
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Using the following relations between r, s and k, x (with assumed to be even): 



-N + r 



1 

1 
2 



X — Xq = < 



s — 



-N + s 

and the following values for Ur and Vs'. 



6k for 1 < r < N/2 

6k for N/2<r<N 

6x for 1 < s < N/2 

6x for N/2 <s<N 



(50) 



Ur 



(51) 



one can show starting from (H71) that fN^x) f{x) as N ^ oo. It's useful to define 



A; 



low 



N -1 TT 



Xlow = Xq — 6x = Xq - [N - 1)^1 TTTT Ax 



2 

iV- 1 



2N 



(52) 



Xhigh = + — ^ — 6x = xo + iN - l)J^Ax 



because then the evaluations of f{k) occur at points, spaced by 6k, the first of which is 
k\ow and the last of which is /chigh; and likewise fN^x) is defined at points spaced by 6x, 
the first of which is xiow and the last of which is Xhigh- 
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